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We present a detailed study of isospin dynamics in peripheral col- 
lisions at Fermi energies. We consider symmetric and mixed colli- 
sions of > Sn isotopes at 35 and bOAMeV to study the isospin 
transport between the different reaction components (residues, gas 
and possibly intermediate mass fragments) and, in particular, the 
charge equilibration in the mixed system. We evaluate the effects 
of drift terms due to asymmetry and density gradients, which are 
directly related to the poorly known value and slope of the sym- 
metry energy below saturation density. We verify the importance 
of an isoscalar momentum dependence of the mean field, which is 
found to influence the isospin transport since it changes the reac- 
tion times. We finally suggest two observables particularly sensitive 
to the isovector part of the nuclear equation-of-state: the correla- 
tion between isospin equilibration and kinetic energy loss for binary 
events, and the isospin content of the produced mid-rapidity frag- 
ments for neck fragmentation events. 
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1 Introduction 



There has been much interest in recent years in the determination of the 
nuclear symmetry energy as a function of density which is of great importance 
for nuclear structure as well as for astrophysical processes. While there has 
been much experimental and theoretical work in this direction in the last years, 
the symmetry energy must still be considered as uncertain. Observables, which 
are sensitive to the isovector density dependence of the Equation-of-State (iso- 
EoS) and testable experimentally are still under investigation [1,2,3,4,5]. 

The symmetry energy E sym appears in the energy density 

e(p, p 3 ) = e(p) + P E sym {p) (p 3 /p) 2 + 0(p 3 /p) 4 + ...., (1) 

expressed in terms of total (p = p p + p n ) and isospin (p 3 = p„ — p p ) densi- 
ties. The symmetry term has a kinetic contribution directly from basic Pauli 
correlations and a potential part from the highly controversial isospin depen- 
dence, as a function of the total nucleon density, of the mean field [4]. Both 
at sub-saturation and supra-saturation densities, predictions based on the ex- 
isting many-body techniques diverge rather widely, see ref. [6]. However, the 
dominant quadratic dependence on the asymmetry parameter (3 = (ps/p) = 
(N — Z)/A in Eq.(l) is well supported by the phenomenology (e.g. the mass 
formula) as well as by all microscopic nuclear many-body calculations [3]. We 
recall that the knowledge of the EoS of asymmetric matter is very important 
at low densities ( e.g. neutron skins, pigmy resonances, nuclear structure at 
the drip lines, neutron distillation in fragmentation, neutron star formation 
and crust) as well as at high densities ( e.g. neutron star mass-radius relation, 
cooling, hybrid structure, transition to a deconfined phase, formation of black 
holes) . 

It is attractive to take advantage of new opportunities in experiments (avail- 
ability of very asymmetric radioactive beams, improved methods of measuring 
event-by-event correlations) and theory (development of reliable microscopic 
transport descriptions for heavy ion collision (HIC)) to obtain results which 
constrain the existing effective interaction models. In this paper we will discuss 
dissipative collisions in the range of Fermi energies, which will yield informa- 
tion on the symmetry term around and below normal density. Here we focus 
our attention on the charge equilibration dynamics in peripheral collisions 
where we expect to see symmetry energy effects on the isospin transport. 
The interesting feature at Fermi energies is the onset of collective flows due 
to compression and expansion of the interacting nuclear matter. The isospin 
transport takes place in regions with density and asymmetry variations and 
thus we expect to have contributions to the isospin current from charge and 
mass drift mechanisms. It was shown [7], that these are determined by the 



2 



value and the density-gradient of the symmetry energy, respectively [7] , which 
are then studied here at subsaturation densities. Fast (pre-equilibrium) nu- 
cleon emission, related to the symmetry energy via the different symmetry 
potentials seen by neutrons and protons, also influences the isospin content of 
the Projectile-like and Target-like residues. Finally in the Fermi energy range 
we have a substantial Intermediate Mass Fragment (IMF) production even in 
non-central collisions and this will also affect the isospin dynamics. We will 
carefully look at all these dynamical effects, trying to select observables which 
are sensitive to the density dependence of the symmetry term. 

Isospin equilibration at intermediate energies and its relation to the symme- 
try energy has, in fact, attracted much attention in recent years in experiment 
as well as in theory [7,8,9,10,11]. We present here a detailed analysis of the 
various iso-transport mechanisms in connection to properties of the effective 
in-medium interactions. We try to disentangle effects due to the general re- 
action dynamics, ruled by isoscalar terms of the interaction, from genuine 
isovector contributions. This will allow us to suggest new observables partic- 
ularly sensitive to the symmetry term of the EoS: 

(i) Correlation between the isospin equilibration of the reaction partners and 
the total kinetic energy loss in binary events; 

(ii) The isospin content of IMF in break-up (ternary) events. 

In Sect. 2 we introduce our transport-theoretical approach to the reaction dy- 
namics and the choice of the effective interactions investigated. In Sect. 3 we 
show results for binary events in Sn + Sn peripheral collisions at 50 AMeV, in 
Sect. 4 we discuss neck fragmentation events, and finally in Sect. 5 we present 
a summary and some perspectives. 



2 Theoretical Description of the Collision Dynamics 

2. 1 The Transport Model 

We perform ab initio collision simulations using the microscopic Stochastic 
Mean Field (SMF) model. It is based on mean field transport theory with cor- 
relations included via hard nucleon-nucleon (NN) collisions and with inclusion 
of stochastic forces acting on the mean phase-space trajectory [4,12,13,14,15]. 
Stochasticity is essential in order to allow the growth of dynamical instabilities 
with fragment production, as well as to obtain physical widths of distributions 
of observables. Moreover it will allow to perform event-by-event correlation 
studies of great importance for the very complex reaction dynamics in this 
energy regime. 
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The transport equation for the phase space distribution function, with the 
Pauli blocking consistently evaluated, is integrated following a representation 
in terms of test particles of finite widths [12,13]. A detailed description of 
the procedure is given in ref. [4]. Our code [16] has been extended by the 
introduction of momentum dependent mean fields (see next subsection), which 
are rather important in this energy range. It has also been possible to improve 
the numerical accuracy while even reducing the computing times [17]. 



A parametrization of free nucleon-nucleon cross sections is used, with isospin, 
energy and angular dependence [18]. Low energy NN collisions, mostly forbid- 
den because of the Pauli blocking, have large cross sections and could induce 
spurious effects in the presence of some numerical inefficiency in the blocking 
procedure, due to the discretization of the phase space. In order to avoid such 
problems a cutoff value a cut = 50mb is used in our calculations. A parallel 
ensemble method is employed in the implementation of the collision term. 



For discussions of isospin dynamics in this energy regime it is essential to have 
a reliable procedure for fragment recognition, i.e. to identify the "gas" (emit- 
ted nucleons and light clusters) and the "liquid" (fragments, residues) phases. 
Fragments are produced via instabilities in the expanding dilute matter and 
as residues of the initial colliding ions (in the case of non-central collisions). 
In our transport simulations the fragment definition is performed by means 
of density cuts ; more precisely, a criterion of p > p cut (with p cut between p/6 
and p/10) for the liquid, and correspondingly for the gas, is used. This pro- 
cedure is applied in an event by event analysis at the "freeze-out", i.e. when 
the resulting fragments are well separated in space and interacting only via 
Coulomb forces. We have checked that such "final" gas/liquid isospin proper- 
ties are not depending on the choice of density cuts and freeze-out times, see 
Sect. 3. 2. Nucleons can be emitted during the early stages of the reaction due 
to hard NN collisions as well as to repulsive potentials, during the fragment 
formation due to the isospin distillation effect, and finally from sequential de- 
cay of excited primary products (fragments/residues). The isospin is used in 
all these dynamical paths both as a tracer of the reaction mechanism, as well 
as an observable of interest with respect to the iso-EoS. 



We follow the reaction dynamics only up to a "freeze-out" time. Fragments 
and residues will still be highly excited and will undergo sequential decays that 
can modify the original isospin information related to the nuclear dynamics. 
However, we will analyze mainly quantities which are expected to change little 
due to the late evaporative emission. "Imbalance" and "Double" Ratios have 
been suggested to that purpose [8,19,20], and will be discussed below. 
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2. 2 Specification of Interactions 



We adopt a generalized form of effective interaction, which can be easily re- 
duced to Skyrme-like forces, with momentum dependent terms also in the 
isovector channel [17,21,22]. The general structure of this isoscalar and isovec- 
tor Momentum Dependent (MD) effective interaction is derived via an asym- 
metric extension of the Gale-Bertsch-DasGupta (GBD) force [23,24,25,26], 
which corresponds to a Yukawian non-locality 
The energy density is parametrized as follows: 

e = e km + e A (A',A")+e B (B',B")+e c (C',C") , (2) 
where Ekin is the usual kinetic energy density and the potential terms are 



e A (A',A") = (A' + A"p 2 )^- 
e B (B',B") = (B' + B"(3 2 ) (j-J g 

e c (C', C") = C'{X NN + 1 PP ) + C"X NP (3) 

The variable (3 = (N — Z)/ (A) defines the isospin content or asymmetry of the 
system, given the number of neutrons (iV), protons (Z), and the total mass 
A = N + Z; g is the nuclear matter density (go is the saturation value). The 
momentum dependence is contained in the X TT > terms, which are integrals of 
the form 

X TT , = J dp dp'f T (f , p) f T > (f, p')g(p, p') , (4) 



with t = P,N, for protons and neutrons. Here f T {^iP) are the nucleon phase 
space distributions for protons and neutrons and the function g(p,p') = 
g[(p — p') 2 ] determines the type of momentum dependence. A Skyrme-like 
momentum dependence is obtained when we use the simple quadratic form 
g{p,p') = {p — p') 2 - A more general momentum dependence, in better agree- 
ment with phenomenological optical potentials, can be introduced by the func- 
tion [23,24,25,26] 



g(p,p') 



i + 



p — p 

A 



(5) 



We remark that this form is particularly suitable for SMF simulations. 
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From the energy density one derives the mean field potentials as U T (f,p) = 
Se/5f T . Thus the above energy density implies a momentum dependent mean 
field interaction. The momentum dependence is isoscalar if the coefficients C 
and C" are identical and also it can have an isovector part, if they are different. 
The isovector momentum dependence implies different effective masses for 
protons and neutrons given as ^ = (1 + j^^f)^ 1 - 

We note that the used form of the Momentum Dependent effective interaction, 
Eqs.(2 - 5), is equivalent to the MDI force introduced in ref.[26]. The choice 
of the parameters is different: having fixed the isoscalar properties of the EoS 
(soft symmetric matter -Kjvm(po) = 215MeV, with nucleon effective mass 
m*/m = 0.67), we are exploring the dynamical effects of a different density 
dependence of the symmetry term, isovector part of the EoS, see Eq.(l). In 
Table I we report the used parametrization. 



Tab. I: Parameters of the Momentum Dependent (MD) Interaction. 



Isoscalar 




Isovector 


Asy- stiff 


Asy-soft 


A' 


-55.626 MeV 


A" 


-41.312 MeV 


59.742 MeV 


B' 


63.013 MeV 


B" 


29.826 MeV 


-71.247 MeV 


a 


-129.41 MeV 


C" 


-1354.78 MeV 


-1354.78 MeV 


a 


1.242 








A 


2.106 fm- 1 









Here we want to test the sensitivity of isospin transport observables to two 
essentially different behaviors of the symmetry energy around saturation: asy- 
soft with a smoothly decreasing behavior below saturation, and asy-stiff with 
instead a rapid decrease to lower densities [1,4]. In Fig.l we show the density 
dependence for these two typical choices. We remark that this is the total 
symmetry energy, with both potential and kinetic contributions. 

At Fermi energies we have substantial density variations during the collision. In 
particular during the expansion phase the system probes more dilute regions. 
Then we expect important effects on the isospin transport due to different 
subsaturation values (larger for asy-soft) and to different slopes (larger for 
asy-stiff) of the symmetry energy. We will see that the transport coefficients 
due to density and isospin gradients (usually called "drift", resp. "diffusion" 
coefficients) will be affected in a different way and the final result will be 
sensitive to the choice of the two E sym shown in Fig.l. 

In the framework of our parametrization of the effective interaction we can 
also easily adjust the isospin momentum dependence to have the same den- 
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Fig. 1. Density dependence of the symmetry energies used in the simulations pre- 
sented here: asy-soft (solid) and asy-stiff (dashed). 

sity dependence of the symmetry energy but opposite neutron/proton effective 
mass splittings. The ordering of the neutron/proton effective masses, which is 
also related to the slope of the Lane potential with energy, is much debated 
in the literature [4,6]. E.g. one could have m* > m* as predicted by the early 
Skyrme forces [27] or the opposite as for the later Skyrme-Lyon parametriza- 
tions [28]. Thus one can, separately from the density dependence, study the 
corresponding dynamical effects on nucleon emissions [17,21,22]. 

The effective mass splitting is increasing with density and it is relevant for the 
dynamics mainly at high momenta. Indeed, it has been shown that the effect on 
nucleon emission and flows can be observed only at high transverse momenta 
[17,21,22]. At the energies investigated here it is, however, a secondary effect, 
as the transport phenomena are mainly driven by the value and the slope of 
the symmetry energy below saturation density. Thus, in order to simplify the 
presentation, we have fixed the isovector momentum dependence. The results 
presented here are obtained with the choice m* > m*. We have actually 
verified that the influence of the different mass splittings on the present results 
on isospin transport is very small. 



2.3 Symmetry Potentials 

Many investigations have given definite evidence that the effective mean fields 
are momentum dependent, i.e. have an important isoscalar momentum de- 
pendence (see e.g. ref. [29]). Only in this way it was possible to explain the 
dependence with energy of the nucleon flow in heavy collisions. We have then 



ASY-SOFT, (3=0.2 ASY-STIFF, (3=0.2 




Fig. 2. Symmetry potentials for neutrons (upper curves) and protons (lower curves) 
in a matter with asymmetry (3 = 0.2 ( 124 Sri-like). Left panel: asy-soft choice; right 
panel: asy-stiff choice; MI interaction: dashed lines, MD interaction: solid {k = /cp) 
and dotted (k = 2kp) lines 

naturally included a momentum dependence of the isoscalar interaction ad- 
justed to these data, as discussed before. However, it is still of interest to 
see the effects of such isoscalar momentum dependence on isospin transport 
phenomena, since this has also been an issue in earlier calculations [11], and 
has not been looked at systematically so far. Thus in the following we will 
also show results with Momentum Independent (MI) interactions along with 
the Momentum Dependent (MD) ones, of course with parameters adjusted 
to the same saturation properties and symmetry terms. Details about the MI 
interaction can be found in ref. [4]. 

We remark that the corresponding symmetry potentials, i.e. the effective field 
variations seen by neutrons and protons in asymmetric matter, can be rather 
different in the two frames, MI vs MD, and this could affect the isospin trans- 
port properties. So this point deserves some more attention. 

Just to have an idea of the expected transport effects we plot in Fig. 2 the 
symmetry potentials for asymmetry matter with (3 = 0.2, the mean asym- 
metry of a l2A Sn isotope. We note that the momentum dependence, in fact, 
modifies the effect of the stiffness of the symmetry term on the nucleon poten- 
tials. E.g. comparing the asy-soft vs. asy-stiff choices for neutrons and for MI 
interactions, we see the general trend of a larger repulsion below saturation 
and a larger attraction above p (and opposite for protons). It is also seen, 
that the effect can be different in the MD case, though depending strongly on 
the nucleon momenta. In this respect we remark that, in particular below the 
saturation density, the differences become appreciable for nucleon momenta 
around 2kp, not much present in the Fermi energy range of interest here. 
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Fig. 3. Density contours of a binary (upper panel) and a ternary (lower panel) event 
in the simulation of 124 5n+ 124 Sn semiperipheral (b = 6/m) collisions at 50AM eV. 
A momentum dependent interaction (MD) has been used. 

At variance a relevant overall effect of the Momentum Dependence is the 
reduction of the interaction times during the collision with an expected clear 
influence on the isospin diffusion. An important consequence is that in order 
to extract information on the density dependence of the symmetry term we 
are forced to select new observables related to the isospin equilibration that 
should be not very sensitive to isoscalar properties of the effective interactions, 
like the main momentum dependence. We will discuss in detail this point in 
the next sections, suggesting very promising possibilities. 



3 General Features of the Isospin Transport in Peripheral Colli- 
sions 



In this work we investigate peripheral collisions between similar systems with 
different isospin, specifically collisions of different combinations of Sn isotopes. 
In Fig. 3 density contours of two typical examples of such events are shown, 
one in the upper panel resulting in two heavy fragments and nucleons or 
light clusters (binary event), and one in the bottom panel producing a third 
fragment of intermediate mass from the neck region between the two heavy 
residues (ternary event, neck fragmentation ). We will discuss binary events in 
Sect. 4, especially with respect to isospin equilibration between the residues, 
and ternary events in Sect. 5 with respect to the isospin content of the neck 
fragments. 
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3.1 Phases of the Reaction Mechanism 

The reaction mechanism can be divided into different phases, which will be 
important for the following discussion: 

(i) Pre-equilibrium phase: the approach and early interaction phase is char- 
acterized by an emission of particles and light clusters (pre-equilibrium 
emission), which contributes to the "gas" phase. Of interest is the asym- 
metry of the gas, which has recently been discussed in detail in central 
collisions in similar systems, in particular, with respect to ratios of pro- 
duced particles and fragments [19,20]. 

(ii) Transport phase: a neck of density below normal values develops between 
the two heavy residues, the evolution of which is driven by the motion of 
the spectators. During this phase isospin is transferred to the neck due 
to the density difference between the neck and the residues; this effect 
has been called isospin migration (isospin drift), which leads to a more 
neutron-rich neck [30,31]. Moreover, in collision systems with different 
asymmetry, isospin is also transported through neck due to the N/Z 
"concentration" difference, leading to an equilibration of the isospin of the 
residues (projectile- and target-like fragments; PLF/TLF). This has been 
called isospin diffusion. In asymmetric systems there is a competition 
of isospin migration and diffusion, which contains information on the 
characteristic time scales of the processes. 

(iii) Fragmentation into identified fragments (primary fragments): as in Fig. 3 
we distinguish two types of events, binary events with two heavy residues 
and the gas, and "ternary" events with an additional fragment of charge 
Z > 3 and with mass between the masses of the residues and the gas, i.e. 
an intermediate mass fragment (IMF). Events with several IMF's are very 
rare at these energies and impact parameters and will not be considered. 

(iv) Evaporation stage: the primary fragments are considerably excited and 
deexcite on the way to the detector, also changing their asymmetry. This 
part is not included in our dynamical model, but could rather be accessed 
with statistical evaporation codes. This is not done here, because we 
try to use observables, which are expected not to change substantially 
in the secondary evaporation. It should be kept in mind, however, that 
the experimental "gas" also contains contributions from the evaporation 
stage. 



3.2 Freeze-out Times and Density Cuts 

In the next sections we will present a detailed study of the isospin structure of 
the reaction products analyzed at the freeze-out times, i.e. of the "primary" 
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Fig. 4. Time evolution of the "gas" asymmetry in the simulation of 124 Sn + 124 Sn 
semiperipheral (b = 6/m) collisions at 50AMeV. Asy-soft: dashed lines. Asy-stiff: 
solid lines. The panels correspond to different density cuts for the "gas" definition, 
see text. A momentum dependent interaction (MD) has been used. 



nuclear systems (nucleons, clusters and fragments) emerging from the reaction 
region. Here we would like to show that the extracted isospin information is 
reliable, i.e. not much dependent on not well denned parameters, like the choice 
of the freeze-out times and of the density cuts needed to separate the liquid 
and gas phases. 

In Fig. 4 we show the time evolution of the "gas" asymmetry averaged over 200 
event simulations of the 124 Sn + 124 Sn semiperipheral (b = 6/m) collisions at 
50AMeV, with the two choices of the symmetry stiffness (MD interactions). 
From the density contour plots of Fig. 3 we see that the freeze-out time should 
be between 150 and 200 fm/c. The flat behaviour of the gas asymmetry after 
150 fm/c clearly indicates that this quantity will not be much sensitive to a 
freeze-out choice in that interval. Of course the same will happen to the liquid 
phase. The calculation has been performed using a po/6 (left panel) and a 
Po/10 (right panel) density cut for the gas phase. We can see that almost 
identical results can be obtained. An interesting point to note is the larger 
final gas asymmetry with the asy-soft choice. This is a combined effect of the 
neutron distillation during the fragment formation (which gives the crossing 
at around 70 fm/c) and of the larger repulsive potential below saturation, see 
Fig.2. 

Tab. II: 124 Sn + 124 Sn collisions: Freeze-out times (fm/c) vs. beam energy, 
centrality and momentum dependence. 
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35 AMeV 


b (fm) 


MD 


MI 


50 AMeV 


b (fm) 


MD 


MI 




6 


180 


200 




6 


140 


150 




8 


160 


175 




8 


120 


125 




10 


120 


125 




10 


100 


100 



Finally in Table II we give our quantitative estimates of the variation of the 
freeze-out times with beam energy, centrality and momentum dependence of 
the used effective interactions. The nuclear interaction times are different for 
MD and MI interactions (shorter for MD) and this property of the isoscalar 
dynamics will certainly affect the isospin diffusion. In order to extract the 
role of the isovector part of the interaction we have to single out observables 
directly related to the interaction times. 



3.3 Isospin Transport Coefficients 



As already pointed out, in heavy ion collisions at Fermi energies isospin trans- 
port is due to the combined effect of asymmetry and density gradients. It can 
be discussed in a compact way by means of the chemical potentials for protons 
and neutrons as a function of density p and isospin asymmetry (5 [7]. In fact, 
the p/n currents can be expressed as 

jp/n = D p p/n V P - D? /n V(3 (6) 



with D p p i n the mass, and D^, the isospin transport coefficients, which are di- 
rectly given by the variation of n,p chemical potentials with respect to density 
and asymmetry, see ref . [7] . In the literature these are often referred to as the 
drift, resp. diffusion coefficients, and we will follow this convention here. Of 
interest here is the differential current of neutrons and protons (iso- vector cur- 
rent) which has a simple relation to the density dependence of the symmetry 
energy. In fact, from the rather general parabolic form of the symmetry term 
in the energy density, we obtain the important relation \i n — p p = Af3E sym and 
derive the drift and diffusion contributions to the isovector current 



Di-D^ApE sym . (7) 

Thus the isospin transport due to density gradients, i.e. isospin migration, 
depends on the slope of the symmetry energy, or the symmetry pressure, 
while the transport due to isospin gradients, i.e. isospin diffusion, depends on 
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the absolute value of the symmetry energy. In peripheral collisions discussed 
here, residues of about normal density are in contact with the neck region of 
density below saturation density. In this region of density a stiff iso-EoS has a 
smaller value but a larger slope compared to a soft iso-EoS. Correspondingly 
we expect opposite effects of these models on the migration and diffusion of 
isospin. 

In our investigation we will see that the asy-soft choice appears to be more 
effective for isospin equilibration. This is the result of a dominant diffusion 
mechanism for peripheral collisions, together with a stronger fast neutron 
emission, which also contributes to the N/Z equilibration. With respect to 
the pre-equilibrium emission we recall that in the asy-soft case neutrons see a 
more repulsive symmetry potential around saturation, see Section 2.3. On the 
other hand, the asy-stiff choice is more effective for isospin migration, which 
will be seen to be important for the isospin content of the neck fragments. 



3.4 The Imbalance Ratio 



We will discuss the asymmetries of the various parts of the reaction system 
(gas, PLF/TLF's, and in the case of ternary events, IMF's ). In particular, we 
study the so-called imbalance ratio (also called Rami or transport ratio [32]), 
which is defined as 

R P,T~ {x H_ xL) > W 



with x eq = \{x H + x L ). Here, x is an isospin sensitive quantity that is to be 
investigated with respect to equilibration. In this work we consider primarily 
the asymmetry (3 = (N — Z)/(N + Z), but also other quantities, such as 
isoscaling coefficients, ratios of production of light fragments, etc, can be of 
interest [5]. The indices H and L refer to the symmetric reaction between the 
heavy (n-rich) and the light (n-poor) systems, while M refers to the mixed 
reaction. P, T denote the rapidity region, in which this quantity is measured, 
in particular the PLF and TLF rapidity regions. Clearly, this ratio is ±1 in 
the projectile and target regions, respectively, for complete transparency, and 
oppositely for complete rebound, while it is zero for complete equilibration. 

In a simple model we can show that the imbalance ratio mainly depends on 
two quantities: the strength of the symmetry energy and the interaction time 
between the two reaction partners. Let us take, for instance, the asymmetry 
(3 of the PLF (or TLF) as the quantity x. In the Fermi energy domain we can 
describe the charge equilibration dynamics as an overdamped dipole oscillation 
mode, see [33,34]. Therefore, as a first order approximation, in the mixed 
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reaction the charge asymmetry parameter will show an exponential relaxation 
towards its complete equilibration value, f3 eq = {(5h + Al)/2 as 

(3™ T = {3 e » + (f3 H ^ L - p eq ) e~ t/T , (9) 

where t is the time elapsed while the reaction partners are interacting (inter- 
action time) and the damping r is mainly connected to the strength of the 
symmetry energy. In fact, as seen in Eq.(7), isospin diffusivity is directly pro- 
portional to the symmetry energy [7,10]. Inserting this expression into Eq.(8), 
one obtains Rp T = ±e - */ r for the PLF and TLF regions, respectively. 

From this simple result one sees that the imbalance ratio does not depend on 
the difference of asymmetries ((3h — (3l) of the systems considered, at least at a 
first order level, which is due essentially to the normalization to the difference 
((3h — (3l) in the definition of R in Eq.(8). Hence the imbalance ratio can be 
considered as a good observable to trace back the strength of the symmetry 
energy from the reaction dynamics. We will show that the effect of different 
iso-EoS choices on the final imbalance ratios are measurable, in particular, 
when correlated to the total energy loss of the dissipative collision, which sets 
the time-scale of the process. 



4 Binary Reactions 

In this section we discuss binary reaction events, and investigate the asym- 
metries of the various reaction components, i.e. of the gas, the residues, and 
the exchanged particles, as well as the imbalance parameter derived from 
these via Eq.(8). We will consider these quantities for the different isovector- 
EoS's, and for momentum- dependent (MD) and momentum-independent (MI) 
parametrizations of the isoscalar part, see Sect. 2.2. The aim of the analysis 
is also to show that effects due to the Momentum-Dependent term of the nu- 
clear interaction are relevant for the reaction mechanism at Fermi energies, 
and consequently for the isospin equilibration. Here we discuss these quanti- 
ties as a function of the incident energy and impact parameter, in the Sect. 4. 3 
we will introduce the more significant correlation to the kinetic energy loss. 

4-1 Asymmetries of Reaction Components 

We first investigate the behaviour of the asymmetries (3 of the reaction com- 
ponents for the different Sn + Sn reactions for the impact parameters b = 6, 8, 
and 10 fm and for two incident energies. We present our results for an incident 
energy of E lab = 50AMeV for the gas phase in Fig. 5, for the residues in Fig.6, 
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Fig. 5. Asymmetries of the gas in Sn + Sn collisions at incident energy of 
E = 50AMeV for MD (left) and MI (right) interactions for a stiff (solid lines) 
and soft (dashed lines) iso-EoS. The symbols differentiate the different reaction sys- 
tems: mixed PLF (full triangle up), mixed TLF (empty triangle down), symmetric 
heavy (full squares), symmetric light (empty, circles ). Thin horizontal lines denote 
the initial asymmetries of the heavy and light system, respectively. 

and for the exchanged particles in Fig. 7. These figures are arranged in an 
analogous way: the results with the momentum dependent interaction (MD) 
are in the left panel, those with the momentum-independent one (MI) on the 
right. Results for the stiff and soft iso-EoS are given by solid and dashed lines, 
respectively. The symbols distinguish the different systems: mixed HL reac- 
tion PLF (triangle up, full), TLF (triangle down, empty), symmetric heavy 
HH (squares, full) and light LL ( circles, empty) systems. Thus, also for the 
gas we have separated the emitted particles according to their origin from the 
projectile (P) or target (T). The asymmetries of the exchanged particles are 
only given for the mixed system. The initial asymmetries of 124 Sn and U2 Sn 
are indicated as horizontal lines. 

We first note the following general properties of these asymmetries, which will 
be valid for all cases discussed in the following. Of course, all these asymmetries 
are connected by charge conservation. In absolute numbers the mass of the 
exchanged particles changes from about 20 to 5 amu, and for the gas from 
about 40 to 20 amu, for b = 6 to 10/m, respectively. 

- The gas (Fig. 5) is generally more n-rich (i.e. asymmetric) than the initial 
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Fig. 6. Asymmetries of the residues in Sn + Sn collisions at incident energy of 
E = 50AMeV for MD (left) and MI (right) interactions. The arrangement of the 
figure and the meaning of the lines and symbols is as in Fig. 5. 

asymmetry of the corresponding nucleus; always for the n-rich nucleus (P) 
but mostly also for the n-poor (less n-rich) system (T). In a nuclear medium 
with neutron excess the interaction is more repulsive for neutrons, such that 
these are emitted preferentially. 

- The asymmetry of the residues (Fig. 6) for the mixed HL system decreases 
for the n-rich (PLF) and increases for the n-poor (TLF) partner with re- 
spect to the initial asymmetries, as expected for isospin equilibration. It is 
interesting to look closer at the residue asymmetry fi Res in the case of the 
HH and LL symmetric collisions. Here we cannot have isospin transport 
and the only variations come from nucleon emissions. The asymmetry de- 
creases for the symmetric systems from the initial asymmetry, less evident 
in the n-poor system, in correspondence with the asymmetry of the gas. 
The change in the asymmetries of the residues of the mixed HL systems is 
much bigger than in the symmetric systems, indicating the importance of 
the transport contributions to the isospin current. This is also supported by 
the dependence on the iso-EoS discussed below. 

- In the mixed HL system neutrons are predominantly exchanged from the 
n-rich to the n-poor system, i.e. the asymmetry of the exchanged particles in 
this direction is much larger than the initial asymmetry of the n-rich system. 
In the other direction the asymmetry is close to the initial asymmetry. 
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Fig. 7. Asymmetries of the exchanged particles in the collision Sn + 112 Sn at 
incident energy of E = 50AMeV for MD (left) and MI (right) interactions, and 
for the stiff (solid lines) and soft (dashed lines) iso-EoS, respectively. PT denote 
particles transferred from the projectile to the target (triangle up, full), and vice 
versa for TP (triangle down, empty). 

- The impact parameter dependence clearly shows that iso-EoS effects are 
more relevant for more dissipative collisions, i.e. for smaller impact parame- 
ters and thus for longer interaction times. This suggests the use of an event 
selection in terms of the total kinetic energy loss in order to enhance the 
sensitivity to the symmetry energy, as discussed in Sect. 4. 3. 

Next we note the differences when using the MD or MI interactions in the 
isoscalar sector. It is seen that the gas is generally more symmetric for a 
MD interaction. The MD interaction has a stronger isospin-blind repulsion, 
which emits protons and neutron equally. Thus the mass of the emitted gas 
is about 5 to 10 units higher in the MD case. This reduces the relative effect 
of the isovector interaction. Correspondingly, the residues are generally more 
symmetric for the MI relative to the MD case since more of the neutron excess 
is emitted into the gas. In addition, in the MD case the reaction proceeds more 
rapidly, because of the greater isoscalar repulsion, thus reducing the interaction 
time and reducing the amount of isospin equilibration. This is seen clearly in 
the mixed systems. 
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With respect to different iso-EoS's we see that the gas is more asymmetric 
for the asy-soft case, which is expected because of the higher sub-saturation 
symmetry energy and thus the larger neutron repulsion. Consequently the 
residues are more symmetric for the soft iso-EoS both for the symmetric and 
the mixed systems. Also, the change of the asymmetries of the residues is much 
larger for the mixed (HL) than for the symmetric systems, indicating that 
isospin equilibration will be mostly sensitive to the iso-EoS via the transport 
contributions of the isospin current. In the exchanged particles we see almost 
no effect for the particles transferred in the direction projectile to target, and 
an interesting difference in the other direction. This is due to the fact that, 
in the direction from projectile to target, the larger isospin diffusion effects 
expected for the asy-soft EoS are counterbalanced by the more neutron-rich 
emission to the gas, which reduces the projectile asymmetry. On the other 
hand, diffusion effects act in opposite direction for the n-poor source, with a 
smaller transferred asymmetry from target to projectile in the asy-soft case, 
also due to the larger neutron emission to the gas. Thus in this direction we 
do expect a sensitivity to the iso-EoS, as seen in Fig. 7. 

However, we have to note that generally the changes due to different isovec- 
tor EoS's are not very large. As we will see below, the imbalance parameter 
considerably increases the sensitivity, since the transport contributions will 
be enhanced. In fact, it was shown in ref. [7] that the transport ratios de- 
pend mostly on (/3pr — Ptp)- i-e. the asymmetry difference of the exchanged 
nucleons, and less on pre-equilibrium emissions. 

Finally in Fig. 8 we compare (for the case of MD interactions) the residue 
asymmetries for incident energies of 50 AMeV (left panel, already shown in 
Fig. 4) and 35 AMeV (right panel). Competing effects are observed, since on 
one hand the emission to the continuum is lower at the lower energy, leading 
to less change of the initial asymmetry, and on the other hand, the interaction 
time is larger, leading to more equilibration in the mixed system. 



4-2 Isospin Imbalance Ratios 

From the asymmetries shown before one obtains the imbalance ratios accord- 
ing to Eq.(8). These measure the amount of isospin equilibration in the mixed 
relative to the symmetric reaction systems, having lower absolute values for 
more equilibration. 

The results for the imbalance ratios are shown in Fig. 9 for the incident energy 
of 50 AMeV (left panel) and 35 AMeV (right panel), for MI/MD interac- 
tions and different iso-EoS choices. Two very clear trends with respect to the 
amount of equilibration appear: the equilibration is larger (R P smaller) both, 
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Fig. 8. Asymmetries of the residues in Sn + Sn collisions at incident energies of 
E = bOAMeV (left) and 35AMeV (right) for MD interactions. The meaning of the 
lines and symbols is as in Fig. 3. 



for lower energies (right panel) and for MI interactions (dashed curves). In 
both cases this can be understood from the fact that the reaction is slower 
and thus the interaction time longer, leading to more equilibration; for the 
lower energy because of the lower speed, and for the MI interaction because of 
the less repulsive isoscalar mean field. In comparing the two iso-EoS's we see 
that the equilibration is larger for the soft iso-EoS, since the higher symmetry 
energy leads to a larger diffusion contribution to the isospin current, as dis- 
cussed above. The first two effects (incident energy and MD interactions) are 
essentially kinematical, depending on the interaction time, while the last one 
depends on the iso-EoS, which is of interest here. Thus, in the next section we 
will propose a scheme to separate the kinematical from the isospin effects. 

The results in Fig. 9 can be compared to other results in the literature. In 
our previous work [7] we have used MI interactions and slightly different 
parametrizations of the iso-EoS, with a somewhat softer symmetry term. The 
results here agree well with the previous ones within these variations. 

On the other hand, there has been extensive discussion of isospin imbalance 
ratios in several papers of B.A. Li and collaborators [11], in which a different 
family of iso-EoS's has been used, which is characterized by an index x (for the 
density dependence of the symmetry energy around saturation of the form p 7 ), 
together with momentum-dependent interactions, and with free and medium- 
modified cross sections. The case x = — 1, corresponding to 7 ~ 1, and with 
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Fig. 9. Imbalance ratios for Sn + Sn collisions for incident energies of 50 (left) and 
35 AMeV (right) as a function of the impact parameter. Signatures of the curves: 
iso-EoS stiff (solid lines), soft (dashed lines); MD interaction (circles), MI interaction 
(squares); projectile rapidity ( full symbols, upper curves ), target rapidity ( open 
symbols, lower curves ). 



free cross sections, should be compared with our stiff MD parametrization 
(7 = 1). However, the value given in ref. [11] for b = 6fm of Rp ~ 0.42 does 
not agree well with our value of Rp ~ 0.25 in Fig. 9. The reason could reside 
in ingredients of the BUU transport code used (point test particles, surface 
effects) leading to a different treatment, relative to the SMF simulations here, 
of the isoscalar part of the nuclear interaction. This will influence, in turn, 
the interaction time between the two reaction partners and the degree of 
dissipation reached in the collision. 

The existing imbalance ratio measurements for Sn isotopes have been analyzed 
with an impact parameter selection performed in terms of charged particle 
multiplicities [8]. A value Rp ~ 0.45 is extracted in the impact parameter 
range of b = 6 — 7.3/m, which is also not far from our asy-stiffMD estimations 
(7 = 1 or slightly stiffer). From Fig. 9 (left panel) we see that in this impact 
parameter range the imbalance ratio Rp is rapidly increasing and so it is very 
important to assess the same centrality selection in the comparison between 
data and simulations. To overcome these problems and isolate isospin effects, 
one could study the imbalance ratio directly as a function of the interaction 
time (or an observable directly related to it), as we will do in the next section. 
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4-3 Correlation with Kinetic Energy Loss 



In subsection 4.1 we have noted, e.g. in Fig.6, that the residues appear more 
symmetric for the momentum-independent interaction. Also, we observe in 
Fig. 8 that at lower energy the longer equilibration time is counterbalanced 
by the smaller pre-equilibrium emission. These trends are more clearly seen 
in Fig.9 for the imbalance ratio, which directly measures the amount of equi- 
libration. In fact the latter is clearly less effective when either the incident 
energy is higher or momentum-dependent interactions are used. As we have 
remarked, both effects decrease the contact time of the two nuclei in the bi- 
nary collision. Thus these observations point to the interaction time as the 
dominating influence on the amount of isospin equilibration, as was already 
discussed in ref. [7]. 

On the other hand, longer interaction times should be correlated with larger 
dissipation and, in particular, with damped PLF and TLF velocities. The 
dissipation, in turn, has been measured, e.g. in deep inelastic collisions, by the 
kinetic energy loss. Thus in this subsection we will investigate the correlation 
between the imbalance parameter and the kinetic energy loss. We define the 
kinetic energy loss per particle as 



where Ap, At, Aplf, Atlf are the masses of the initial projectile and target, 
and of the final projectile-like and target-like fragments, respectively. Here 
E cm is the initial energy, per nucleon, available in the cm system. Ekm and 
Ep° t ul are the final kinetic energies of the fragments in the cm system, and the 
Coulomb potential energy, respectively. We will use the relative energy loss 
Eioss/Ecm as a measure of PLF/TLF velocities and thus of the interaction 
time. The study of isospin equilibration as a function of the heavy residue 
excitation energy (related to the kinetic energy loss ) was also suggested in 



The relative energy loss is shown in Fig. 10 (left panel) as a function of impact 
parameter for the two incident energies and MD, resp. MI interactions. It 
is seen that this quantity exhibits, generally, a very good correlation with 
impact parameter, except for the large impact parameters, where eventually 
all curves will converge to zero. As expected the energy loss at a given impact 
parameter is larger for lower energies and for MI interactions, because of the 
longer interaction times. 

However, the energy loss consists not only in internal excitation of the residues 
but also in kinetic energy of the particles emitted to the gas. The mass of the 
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ref. [9]. 
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Fig. 10. (left panel) Relative energy loss as a function of impact parameter; and 
(right panel) mass of particles emitted into the g function of relative energy 

loss. Curves are shown for MD interactions (solid lines and circles) and MI interac- 
tions (dashed lines and squares), and for 35 MeV (open symbols) and 50 MeV (full 
symbols) incident energy, respectively. 

gas is shown in the right hand panel of Fig. 10 as a function of the relative en- 
ergy loss. One sees that the mass loss (and correspondingly the energy loss) to 
the gas is relatively independent of incident energy for MD interactions (black 
curves), while there is a larger difference for the MI interactions (grey/red 
curves). Since the energy lost into the gas is not necessarily closely related to 
the interaction time, this represents a correction to the correlation between 
interaction time and energy loss, which is discussed below. 

In Fig. 11 we present the imbalance ratio as a function of the relative energy 
loss. In the upper panel we separate the results for MD and MI interactions 
and stiff and soft iso-EoS, but in each case the results for 35 and 50 MeV 
are collected together and connected by lines. In contrast to Fig. 9 we now 
see that the points for the different incident energies approximately fall on 
one line, however with considerable scatter around it. The deviations from a 
smooth behaviour are larger for the MI interactions (diamonds and triangles). 
This can be traced back to the observation, made above, that in this case 
the pre-equilibrium emission, i.e. the energy loss to the gas, is rather different 
for the two beam energies, while this is not so much the case for the MD 
interaction. Considering that the pre-equilibrium emission is not a measure of 
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Fig. 11. Imbalance ratio as a function of relative energy loss. Upper panel: Separately 
for stiff (solid) and soft (dashed) iso-EoS, and for MD (circles and squares) and MI 
(diamonds and triangles) interactions, in the projectile region (full symbols) and 
the target region (open symbols). Lower panel: Quadratic fit to all points for the 
stiff (solid), resp. soft (dashed) iso-EoS. 

the interaction time, this introduces a correction to the values obtained for 35 
MeV relative to those at 50 MeV, which tends to make the dependence more 
smooth. We do not attempt, however, to make this correction quantitative, 
because this would amount to separating the gas particles into pre-equilibrium 
and evaporation components. 

It is seen in Fig. 11 that the curves for the asy-soft EoS (dashed) are generally 
lower in the projectile region (and oppositely for the target region), i.e. show 
more equilibration, that those for the asy-stiff EoS. In order to emphasize 
this trend we have, in the lower panel of the figure, collected together all the 
values for the stiff (circles) and the soft (squares) iso-EoS, and fitted them by 
a quadratic curve. It is seen that this fit gives a good representation of the 
trend of the results. 

The difference between the curves for the stiff and soft iso-EOS in the lower 
panel then isolates the influence of the iso-EoS from kinematical effects asso- 
ciated with the interaction time. It is seen, that there is a systematic effect of 
the symmetry energy of the order of about 20 percent, which should be mea- 
surable. The correlation suggested in Fig. 11 should represent a general feature 
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of isospin diffusion, and it would be of great interest to verify experimentally. 
As discussed above, this kind of analysis would help also in the comparison 
with the results of other theoretical models. 

On the other hand, it is also seen that the effect is not extremely large, in view 
also of the substantial scatter of the points around the curves. Thus it is still 
of interest to investigate other observables, which might show more sensitivity 
to the iso-EoS. This will be done in the next section with respect to ternary 
events. 



5 Ternary Reactions 

5. 1 Characterization of Intermediate Mass Fragments 

As discussed above, and seen in the example of Fig. 3, there are peripheral 
events in which a third intermediate mass fragment (IMF) appears. The term 
"neck fragmentation" has been coined for such events, see [35] and refs. therein. 
Related isospin transport studies have been performed theoretically for sym- 
metric reactions [7] and theoretically and experimentally for asymmetric col- 
lisions [31,36]. The asymmetry and the kinematical properties of the IMF 
were shown to carry additional information on the symmetry energy. Here we 
investigate this information more in detail in the system under study. 

In such reactions the asymmetry of the IMF is of particular interest. In sym- 
metric systems it is a result of isospin migration, due to density gradients 
between the residues and the dilute neck. In mixed systems this effect is in 
competition with the isospin diffusion discussed in the last section due to con- 
centration gradients. While the study of this competition may yield further 
information on the time scales of these different processes, the more pure case 
will be the study of the symmetric systems. 

We start with characterizing the origin of the IMF's. As an example in Fig. 12 
we display various properties of IMF's from the mixed reaction 124 Sn+ 112 Sn at 
50 AMeV and an impact parameter of b = Qfm for the stiff (solid histogram, 
black symbols) and soft (dotted, grey/red ) iso-EoS. In the upper left panel 
the charge distribution is displayed, exhibiting a broad distribution around 
charge Z ks 8, which is not very different for different iso-EoS's. The upper 
right panel shows the times at which the third fragments are identified in 
our method. It is seen that they appear rather late and in a rather short time 
span. The lower left panel displays the asymmetry of the IMF's as a function of 
their charge (later we will only show the average asymmetry). This asymmetry 
distribution is rather flat; it is, however, substantially different depending on 
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Fig. 12. Properties of Intermediate Mass Fragments (IMF) in the reaction 
124 Sn+ 112 Sn at 50 AMeV and b = 6fm for the stiff (solid lines and black symbols) 
and soft (dotted lines, grey/red symbols) iso-EoS: (upper left) charge distribution; 
(upper right) time of appearance of IMF; (lower left) asymmetry distribution of 
IMF's; (lower right) velocity correlation of IMF with respect to PLF (r) and TLF 
(rl), where r is the ratio of the relative IMF-residue velocity to the Viola systematics 
velocity. 

the iso-EoS [30]. Finally in the lower right panel we show velocity correlations 
between the fragments and the residues, as it was suggested in ref. [31] (so- 
called Wilczynski-2 plot). The quantity r is the ratio of the relative velocity of 
the IMF and the PLF to the Viola systematics velocity, which corresponds to 
fission-like events, and rl is the same quantity with respect to the TLF. A value 
of 1 signifies a late, statistical origin of the IMF, while values substantially 
different from one indicate an IMF of dynamical origin, which is clearly the 
case here. Altogether the information in this figure then suggests that the 
IMF's come from a uniform source, which can be identified with the neck 
region. 



5.2 Isospin Dynamics with Fragment Production 



We have looked at the asymmetries of the various reactions parts (gas, residues 
and exchanged particles, like in Figs. 5, 6 and 7) also for the ternary events. 
It could be expected that these are different from those of binary events, 
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Fig. 13. Asymmetries of IMF's in ternary Sn+Sn reactions at 50 AMeV as a function 
of impact parameter for MD (left panel) and MI (right panel) interactions: symmet- 
ric 124+124 (squares), symmetric 112+112 (diamonds), mixed 124+112 (circles); 
stiff iso-EoS (solid lines), soft iso-EoS (dashed lines). Horizontal thin lines denote 
the asymmetries of 124 Sn, mixed reaction and 112 5n, respectively. 

because some of the isospin is carried away by the IMF. However, we find that 
the differences are very small and thus we do not show these figures separately. 
The relatively small mass of the IMF's thus does not substantially influence 
the asymmetry of the main reaction products. 

A quantity which is very sensitive to the iso-EoS is the asymmetry of the 
IMF. It is shown in Fig. 13 at 50 AMeV for the MD (left panel) and MI (right 
panel) interactions, separately for the symmetric and the mixed reactions. The 
asymmetry of the IMF is larger, i.e. the IMF is more n-rich, for the stiff relative 
to the soft iso-EoS, since the former exhibits a larger isospin migration due to 
the larger slope of the symmetry energy below saturation. This is clearly the 
case for the symmetric reactions, but it is also true for the mixed reactions, 
where there is a competition with isospin diffusion, which depends on the value 
of the symmetry energy and it is larger for the soft iso-EoS. Our result then 
shows that the isospin migration is the dominating effect for the asymmetry 
of the neck fragments. 

One may now think to investigate the imbalance ratio for the asymmetry of the 
IMF, i.e. the ratio in Eq.(8) where the quantity x represents the corresponding 
average asymmetry of the IMF's independent of their rapidity. As was already 
seen in ref. [30], this quantity is always found to be close to zero. This means 
that the asymmetry of the IMF in the mixed reaction is the average of that in 
the symmetric reactions. This is expected around mid-rapidity for symmetry 
reasons, although the IMF-asymmetries for all the systems, HH, LL and 
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HL, are different from the initial values. This is then approximately still true 
for all the neck fragments. In fact we note that the IMF asymmetries for all 
systems are changing in the same direction, increasing with the stiffness of the 
symmetry term. This is a non-trivial condition for the almost zero imbalance 
ratio of the neck fragments. 

However, a very striking feature seen in Fig. 13 is the fact that the difference 
between the IMF-asymmetry for the stiff and the soft iso-EoS's is large, and 
in addition, that it is much larger in the case of MD interactions relative to 
MI interactions (difference between solid and dashed lines in the two panels, 
for each colliding system). The large sensitivity of the IMF asymmetry to 
the iso-EOS was already noted earlier [30] for MI interactions, showing the 
effectiveness of isospin migration for a stiff iso-EoS. 

An interesting point which deserves a deeper interpretation is that the iso-EoS 
effect appears even more prominent for MD interactions. We see that (3jmf 
increases for the asy-stiff case (solid lines) and decreases for the asy-soft one 
(dashed lines) for MD relative to MI interactions. Therefore the differences 
in the predicted IMF-asymmetries for different iso-EoS's are systematically 
enhanced. This can be traced back to the following effect of MD interactions: 
the reaction is faster and so we obtain a mid-rapidity neck fragment only in 
regions with larger projectile-target overlaps. Thus, for a given impact param- 
eter the IMF is produced in a more compact configuration, i.e. with a larger 
interface with the spectators [37] near normal density, where the asy-stiff EoS 
is effective with respect to neutron migration to the neck, while the asy-soft 
is not, due to the low value of the derivative of the symmetry energy close to 
normal density (see Fig.l). In conclusion, with MD interactions we see a more 
prominent iso-EoS effect with a strong n-enrichment of the neck fragments 
only for the asy-stiff case. 



5.3 Asymmetry Ratios 

It is of interest to construct an observable, which quantifies these effects and 
which does not depend sensitively on measuring the absolute asymmetries, 
which are changed by the secondary evaporation. Considering that the asym- 
metries of the other reaction products (residues and the gas) is not equally 
sensitive to the iso-EoS (see Figs. 5 and 6 ; which are for binary events, but, as 
stated above, do not change noticeably for ternary events), it is attractive to 
look for ratios of asymmetries. In Fig. 14 we show the ratio of the asymmetries 
of the IMF's to those of the residues for the symmetric Sn + Sn reactions for 
stiff (solid) and soft (dashed) iso-EoS, and for MD (circles) and MI (squares) 
interactions. The results correspond to b = 6fm semiperipheral events, plotted 
here as a function of the initial isospin asymmetry of the colliding system. 
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Fig. 14. Ratio of asymmetries of IMF to residues for symmetric Sn + Sn reactions 
as a function of the initial isospin asymmetry for semiperipheral events (b = 6 fin) 
at 50AMeV beam energy. Results are shown for MD (circles) and MI (squares) 
interactions, and for stiff (solid lines) and soft (dashed lines) iso-EoS's. 

The ratio between the asymmetry of IMF's and residues can be estimated on 
the basis of simple energy balance considerations. In fact, isospin migration 
is due to the fact that the neck region has lower density with respect to the 
residues and the symmetry energy is decreasing with density. Starting from a 
residue of mass A res and a neck of mass Ajmf of different density but uniform 
asymmetry /3, we assume that the mass A participating in the isospin exchange 
is approximately equal to the mass of the neck, while it is small relative to the 
mass of the residue. This will lead to the asymmetry (j3 + Aj3) of the neck, and 
to a total asymmetry (3 res = [(5{A res - A) + (/3 - A(3)A]/A res = (3 - A/3A/A res 
of the residue, with A/3 to be determined by minimization of the symmetry 
energy. The corresponding variation of the symmetry energy is equal (apart 
from a constant) to: 

AE sym = A res E sym (p R )(P - A(3A/A res ) 2 + AE sym { Pl )((3 + A(3) 2 , (11) 
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where pr and pi are the densities of the residue and neck regions, respectively. 
The minimum of the variation of AE sym yields 
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From this simple argument the ratio between the IMF and residue asymmetries 
should depend only on symmetry energy properties and, in particular, on the 
difference of the symmetry energy between the residue and the neck regions, 
as appropriate for isospin migration. It should also be larger than one, more 
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Fig. 15. Ratio of asymmetries of IMF to gas for symmetric Sn + Sn reactions as 
a function of the initial isospin asymmetry. Semiperipheral events (b = 6fm) at 
50AMeV beam energy. Signatures are as in Fig. 12 

so for the asy-stiff than for the asy-soft EoS. 

It is seen indeed in Fig. 12, that this ratio of IMF over residues asymmetry is 
nicely dependent on the iso-EoS only and not on the initial asymmetry. This 
is a clear feature of the dominant isospin migration mechanism which depends 
only on the density gradients. The ratio is also mostly larger than unity. In 
both cases, MI and MD interactions, the asy-stiff 'EoS is more effective in the n- 
enrichment of the dilute region where the IMF's are produced. Moreover, since 
for the asy-stiff 'EoS, the derivative of the symmetry energy with respect to the 
density is constant, isospin migration effects should be essentially the same in 
MI and MD interactions, as roughly confirmed by Fig. 12. On the other hand, 
isospin effects appear quite reduced for the asy-soft EoS, especially for the 
MD interaction for the reasons already mentioned (reduced slope and larger 
neutron emission for the asy-soft EoS for the more compact configurations 
of the MD dynamics). These behaviours leads to a rather large sensitivity to 
the iso-EOS for the more realistic MD interactions. We note that the iso-EoS 
effect amounts to more than 30%. Since we expect the secondary evaporation 
to affect in a similar way the asymmetry of IMF's and residues, we suggest 
this observable as very sensitive probe to the isovector part of the EoS. 

A similar investigation can be performed with the ratio of the asymmetries of 
the IMF to those of the gas. The corresponding results are shown in Fig. 15. The 
mechanism of the gas emission is not simply related to the isospin migration 
and now the ratios are seen to be more dependent of the initial asymmetry. 
The decreasing trend seen from 112 Sn to 12A Sn collisions is a consequence of 
the larger asymmetry of the gas for the more isospin asymmetric systems (see 
Fig. 5). Impressive is the large value of this ratio for the 112 Sn collisions and 
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MD interaction and asy-stiff EoS, which nicely confirms the interpretation 
of a larger proton enrichment and thus smaller asymmetry of the gas phase. 
Also for this observable a larger sensitivity to the iso-EoS is observed in the 
MD case. This reflects the large difference in the predicted IMF asymmetry 
between the two iso-EoS's, as discussed above. 

5.4 Double Ratios 

In order to further reduce the secondary decay effects it has recently been 
proposed to look at double ratios of isospin-dependent quantities, n-rich vs. 
n-poor systems, e.g. in central collisions for the asymmetries of nucleon [19] 
or IMF [20] emissions. From Figs. 14/15 this double ratio should always be 
around unity for Pimf/ fines, n ot much depending on the initial asymmetry, 
while it appears of interest for the quantity (3jmfI (3 gas. From the same class 
of semiperipheral events we can extract Tab. Ill for the Double Ratios: 

(3 IMF /f3gas( 1M Sn + 124 Sn) 
/3 IMF /Pgas( 112 Sn + 112 Sn) 

We note the almost 20% Iso-EoS dependence for the more realistic MD in- 
teractions. 

Tab. Ill: Pimf / '13 gas Double Ratio. 



Interaction 


MD 


MI 


Asy — soft 


0.660 


0.687 


Asy - stiff 


0.572 


0.628 



6 Summary and Perspectives 

In this work we have systematically studied isospin transport in heavy ion 
collisions in the Fermi energy domain for symmetric and asymmetric combina- 
tions of Sn isotopes. We have focused on semi-peripheral dissipative collisions 
to study both isospin equilibration between the residues and fragmentation 
of the neck between the residues, i.e. to study the isovector part of the EoS 
for densities below saturation density. The purpose of this investigation has 
been to analyze in detail the mechanisms of isospin dynamics and to iden- 
tify observables, which are sensitive to the still controversial iso-EoS. We have 
employed two typically different iso-EoS's (asy-soft and asy-stiff), where one 
expects to see characteristic signatures since the value and the density-slope of 
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the symmetry energy determine the effects of isospin migration (due to density 
differences) and isospin diffusion (due to isospin concentration differences). 

We have clearly shown that the isospin transport is also dependent on the over- 
all reaction dynamics, mostly ruled by the isoscalar properties of the effective 
interactions. Thus, it is important to select observables able to disentangle 
between isoscalar and isovector contributions. 

In this respect we have thoroughly analyzed the effects of the well established 
(isoscalar) momentum dependence of the mean field on isospin transport phe- 
nomena. The momentum dependence of the effective forces can affect the 
relationship between isospin transport observables and the symmetry energy 
since the reaction dynamics is modified. In particular, the interaction times 
are shorter and the composite nuclear systems break starting from more com- 
pact configurations. It is then difficult to distinguish effects of the iso-EoS and 
of the momentum-dependence in the usual analysis as a function of centrality. 

For binary dissipative reactions we therefore suggest a study of the imbalance 
ratios as a function of the relative energy loss, which is in turn well correlated 
to the interaction time. We thus obtain rather well defined universal curves, 
including results from different MD interactions at different beam energies, 
which present, however, a clear dependence on the iso-EoS of the order of 
20%. This kind of analysis should be useful and illuminating in the compar- 
ison with experimental data, as well as in testing the predictions of different 
theoretical models. In fact, it allows to disentangle isospin effects from other 
model ingredients that may influence the degree of dissipation reached in the 
system. 

For ternary events, where an Intermediate Mass Fragment (IMF) is formed in 
the mid-rapidity region (neck fragmentation), the IMF neutron enrichment is 
very sensitive to the iso-EoS. We have isolated the isospin migration as the 
dominant mechanism in this process in collisions of the symmetric systems 
( 124 Sn+ 124 Sn vs. 112 Sn+ 112 Sn). The iso-EoS dependence is clearly emerging, 
being stronger for an asy-stiff choice with a larger slope of the symmetry 
energy. We have also found that the effect is larger for the more realistic MD 
dynamics, due to different density regions which are probed. 

All the properties of the reaction products, i.e. residues, gas and IMF's, are 
evaluated at the freeze-out time and thus they can be modified by the sequen- 
tial decay of the excited products. We therefore propose the investigation of 
"relative" isospin contents, i.e. of ratios of the asymmetry (3 = (N — Z)/A 
of the various products, which are expected to be less affected by secondary 
emissions. A very promising quantity seems to be the ratio (3imf/ fines at mid- 
centrality, which directly measures the isospin migration. In fact, it is found 
to be very sensitive to the iso-EoS, up to above the 30% level for the more 
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realistic MD interactions. 

In conclusion, we suggest new isospin sensitive observables to study in dissi- 
pative HIC at Fermi energies in order to further constrain value and slope of 
the symmetry term at subsaturation densities. 

It should also be noticed that the ratios considered here (Imbalance Ratios, 
Pimf / fines) do not depend, at a first level approximation, on the initial asym- 
metry of the systems considered, but only on the properties (value and slope) 
of the symmetry energy. Hence valuable information on these properties can 
already be obtained with the available neutron-rich beams. However, the use 
of radioactive beams in this intermediate energy range would be certainly im- 
portant in order to extend the systematics of results and to better constrain 
the symmetry energy behaviour. 
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